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We formulate a semiclassical approach to study the dynamics of coherence loss and revival in a 
Bose-Josephson dimer. The phase-space structure of the bi-modal system in the Rabi, Josephson, 
and Fock interaction regimes, is reviewed and the prescription for its WKB quantization is specified. 
The local density of states (LDOS) is then deduced for any given preparation from its semiclassical 
projection onto the WKB eigenstates. The LDOS and the non-linear variation of its level-spacing 
are employed to construct the time evolution of the initial preparation and study the temporal 
fluctuations of interferometric fringe visibility. The qualitative behavior and characteristic timescales 
of these fiuctuations are set by the pertinent participation number, quantifying the spectral content 
of the preparation. We employ this methodology to study the Josephson-regime coherence dynamics 
of several initial state preparations, including a Twin-Fock state and three different coherent states 
that we denote as 'Zero', 'Pi', and 'Edge' (the latter two are both on-separatrix preparations, while 
the Zero is the standard ground sate preparation). We find a remarkable agreement between the 
semiclassical predictions and numerical simulations of the full quantum dynamics. Consequently, a 
characteristic distinct behavior is implied for each of the different preparations. 



I. INTRODUCTION 

Atom interfcrometry [l|, Q with Bosc-condensed atoms offers the possibility of constructing compact and highly 
precise measurement tools. Recent experiments [sl-flU demonstrate that bi-modal Bose-Einstein condensates (BECs) 
have the necessary phase-coherence and controllability of coupling and interaction parameters, to operate atom inter- 
ferometers at the best sensitivity allowed by quantum mechanics. 

A typical atom interferometer follows the Mach-Zehnder scheme consisting of a preparation stage, in which the 
bimodal input state is mixed by a (usually 50:50) beam-splitter, a waiting time (t) during which the system evolves 
and the two modes acquire a relative phase difference (p, and a measurement stage where the two condensates are either 
released and allowed to interfere or are mixed again by a second beamsplitter. In the former case, the accumulated 
relative phase (p correlates with the location of interference fringes, whereas in the latter it is reflected in the final 
atom number difference. 

The single-particle coherence of the split BEC is characterized by the many-realizations fringe-visibility function 

(7^2^ (t). The useful timescale of an interferometric measurement is set by decoherence. For trapped BEC interfer- 
ometers, an important source for the loss of single-particle coherence, is the phase-diffusion induced by nonlinear 
inter-particle interactions 0, [l6l - l22j . For uncoupled condensates starting from a coherent preparation, this process 

amounts to a simple Gaussian decay of g[^2i^) on a time scale t (C/a/ZV)^^, where U is the interaction strength 
and N is the number of particles. However, most current atom-interferometer setups with trapped BECs operate in 
the Josephson interaction regime, where coupling while small, still affects the dynamics and generates richer coherence 
evolution which includes oscillations, fluctuations, and recurrences of g[^2i^) "^^^ ^o quantum collapse and revival. 

Our objective here is to characterize these fluctuations and explore their dependence on the initial preparation 
j23l - l26j . For this purpose we employ a semiclassical picture of the quantum two-mode Bose-Hubbard model normally 
used to describe BEC atom interferometers. In Section II We introduce the model Hamiltonian for N boson in a dimer 
system. In section III we provide the prescription for the WKB quantization of the associated spherical phase space. 
Several experimentally viable preparation are introduced in Section IV, and their local density of states (LDOS) 
is characterized using a participation number Af , with marked differences in the dependence of M on the particle 
number N, and on the interaction U. Consequently we are able to analyze the observed temporal fluctuations in 
Section V and determine their long time average and their characteristic variance. We show that the RMS of the 
fluctuations scales differently with TV, depending on the nature of the prepared state. 
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II. MODEL HAMILTONIAN AND PHASE-SPACE STRUCTURE 



Matter-wave interferometers of interacting atoms are often described by the tight-binding Bose-Hubbard Hamiho- 
nian (BHH) pol - l22l . [27l - l29j . Here we refer to particles in a two-site (bi- modal) system, also known as a 'dimer': 



i=l,2 



Sifi, + ^fitini - 1) 



- y(a2ai + a,\a2) , 



where K is the hopping amplitude, U is the interaction, and £i are the on-site energies. Defining 

J. = ^(ni - 722) 



1. 

2K 



n. 



J+ = 0102, 



the Hamiltonian ([T]) is formally the same as that of a spin j = N/2 system. 



n 



-£Jz + UJt - KJ.^ 



const 



(1) 

(2) 
(3) 

(4) 



where £ = £2 — £1 is the bias in the on-site potentials. Rescaling time as i — Kt it is clear that the characteristic 
dimensionless parameters which determine both stationary and dynamic properties, are 



NU/K 
£/K 



(5) 
(6) 



The two-site BHH can be regarded as the quantized version of the top Hamiltonian, whose spherical phase space is 
described by the conjugate non-canonical coordinates (9, ip) that are defined through 



■h = [(j + l)jf cos(0), 
Jx = [(j + l)jf/' sin(0)cos(^). 
The hamiltonian (|4]) is thus transformed into the top Hamiltonian, 



NK 



ecosi 



sm t> cos ip 



(7) 
(8) 



(9) 



which is closely related to the Josephson Hamiltonian. The cylindrical phase space of the latter is described by the 
conjugate canonical coordinates n = cos9 and f. In the absence of bias, near the equator, it is just the pendulum 
Hamiltonian: 



H{n,(p) « NK 



un — cos ip 



(10) 



In Fig[T]we draw the constant energy contours T-L{9, ip) = const of the top Hamiltonian © for u = 10 and e — 0. Gen- 
erally, the qualitative features of the phase-space structure [itI . [30| , change drastically with the interaction strength. 
For u > 1 a separatrix appears provided jej < Ec, where 



^2/3 - 1 



3/2 



(11) 



This separatrix divides the spherical phase space into "sea" and two "islands" as shown in FigUl In what follows, 
we focus on the case of zero bias. Accordingly we distinguish between 3 regimes depending on the strength of the 
interaction 1311: 



Rabi regime: m < 1 , 
Josephson regime: 1 < w < A^^ 
Fock regime: u > N'^ . 



(12) 
(13) 
(14) 



In the Rabi regime the separatrix disappears and the entire phase-space consists of the nearly linear "sea". By 
contrast, in the Fock regimes the "sea" has area less than 1/N, and therefore it cannot accommodate quantum states. 
Thus in the Fock regime phase space is composed entirely of two nonlinear components: the "islands" occupy the 
entire upper and lower hemispheres. Our main interest below is in the intermediate Josephson regime where linear 
and non-linear regions coexist and the dynamics is least trivial. This is the regime of interest to most current BEG 
interferometers and luckily, precisely where scmiclassical methods are expected to be most effective. 
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FIG. 1: (Color online) Contour lines for u > 2. Sea levels are colored blue, Island levels are colored green, and the Separatrix is 
colored red (left panel). Energy spectrum for A'' = 20 and u = 10. WKB energies (red x) are compared with exact eigenvalues 
(blue +). Dashed lines indicate slopes ujj for low energies, lj^ for near-separatrix energies, and a;+ for high energies (right 
panel) . 



III. WKB QUANTIZATION AND CHARACTERISTIC FREQUENCIES/LEVEL-SPACINGS 



Due to the simplicity of the two mode BHH ([3]) , it is possible to carry out its semiclassical quantization analyti cally 
[s^ - fs^ l and acquire great insight on the ensuing dynamics of the corresponding Wigner distribution §^ [36| . 
We begin by accurately determining the quantum energy levels in the Josephson regime, using the WKB prescription. 
Having a phase-space area of Att spherical angle supporting iV + 1 quantum states, the Planck cell area is. 



h = Planck cell area in steradians = 



47r 
iV+1 ' 



and the WKB quantization condition thus reads. 



1 

2^^ 



j/ = 0, 1,2,3,, 



(15) 



(16) 



where A{E) is the phase space area which is enclosed by the energy contour. Note that it does not matter which area, 
of which "side" of the contour, is selected. Away from the separatrix the levels spacing equals approximately to the 
classical oscillation frequency: 



dE 
dv 



A'{E) 



(17) 

(18) 
(19) 
(20) 

(21) 

Comparing Eq. (PT|) to Eq. ([T51) and Eq. (pil)) we realize that only in the vicinity of the separatrix does the number 
of particles [N) become an essential parameter in the spectral analysis of the dynamics at fixed u. 



In particular in the absence of interaction this is the Rabi frequency 

UJK ~ K 

For strong interaction, in the bottom of the sea, it is the Josephson frequency 

ujj w VNUK = \/uu}K 

while in the top of the islands it is 

uj-^ « NU = u ujk 
Finally, in the vicinity of the separatrix a more careful analysis is required, leading to the h dependent result 





log 




-1 




m 
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IV. THE PREPARATIONS 



Current experiments in matter- wave interferometry enable the preparation of nearly coherent states of the SU (2) 
algebra, by fast splitting of the condensate [loj . Alternatively, number-squeezed states, approaching relative- number 



Fock states for large separation, can be prepared by slow, adiabatic splitting [12|. So far, the main focus of study 
in the Josephson regime, has been on coherent population dynamics, contrasting coherent preparations located near 
the bottom of the linear sea which exhibit Josephson oscillations around the ground state [g, [37h39| | , with coherent 
preparations located near the top of the nonlinear islands, which result in self-trapped phase-oscillations around the 
'poles' 0, |4^. Our focus here, is on the interesting effects incurred in the fringe-visibility dynamics of coherent 
preparations located on the separatrix and contrasting them with the more common ground-state, north-pole, and 
number-squeezed preparations. 



A. Wigner distributions 



In order to gain semiclassical insight, it is convenient to represent each eigenstates \Eu) by a proper spin Wigner 
function (ssl. [Sa]. which is a quasi-distribution that dwells on the spherical phase space. In this representation these 
eigenstates corresponds to strips along the contour lines oiT-L. In the same representation Coherent states \Otf) are 
like a minimal Gaussian wavepackets, while Fock states \n) are like equi-latitude annulus. Note that the coherent 
state 9=Q is also a Fock state with all the particles occupying one site, while \n = 0) is the Twin Fock state with 
equal number of particles in both sites. 




FIG. 2: (Color online) An illustration of the NorthPole (^=0) coherent state preparation (left), the TwinFock (n = 0) prepa- 
ration (middle), and of Pi ("tt"). Zero ("0") and Edge ("e") preparations (right) using Wigner plots on a sphere. The left and 
middle panels are a 3D plots, while the right panel is a Mercator projection of the sphere using {(p, n) coordinates. 



In Fig.[5] we plot the Wigner functions corresponding to the five preparations under study. These include the 
NorthPole self-trapped state, the TwinFock state, two equal-population coherent states that we call Zero {f—Q) and 
Pi ((/j=7r), and a third coherent state preparation that we call Edge. The two latter states (Pi and Edge) are both 
on-separatrix preparations. Note that in the Zero state all the particles occupy the symmetric orbital, while in the 
highly excited Pi state all the particles occupy the antisymmetric orbital. 

Some of these preparations were experimentally realized and studied @. While in the experimental work the 
emphasis was on contrasting Josephson- Oscillation and Self- Trapping, with regard to near-Zero and near-NorthPole 
preparations, here our main interest is in contrasting the Zero preparation with the on-separatrix preparations Pi and 
Edge. 



B. Local Density of States 



In order to analyze the dynamics ensuing from the said preparations, we expand the initial state as a superposition 
of the eigenstates \Eu). The probability of the i/th eigenstate in the superposition is denoted P(£'i/), and is known as 
the local density of states (LDOS) with respect to the pertinent preparation. The LDOS of the various preparations 
is illustrated in Fig.[3] and the line shape can be determined analytically via a semiclassical calculation [2^ [26| . 
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FIG. 3: (Color online) The LDOS oi N = 500 bosons with u — A, for TwinFock, Zero, Pi, and Edge preparations (left to 
right). The horizontal axes are E — and uj/loj. The lines in the LDOS figures are based on a semiclassical analysis, while 
the circles are from the exact quantum calculation. Note the outstanding difference between the spectral support of Zero and 
Pi preparations compared with continuous-like support in the case of Edge and Fock preparations. 



Schematically the results can be summarized as follows: 

TwinFock 

^ ' Zero 
V{E) 



1 - 



2E 



nE) 



Pi 



Edge 



I 

K 

exp 



NK 

E -E_' 



-I -1/2 



NU 
E-E.. 
NU 

1 fE-E^ 



N \ ojj 



(22) 

(23) 
(24) 

(25) 



where I and K are Bessel functions. It is important to observe that the classical energy scales are NK and N'^U. 
Accordingly only the line shape of the TwinFock LDOS has a purely classical interpretation. In contrast to that, the 
width of the coherent preparations is determined by the quantum uncertainty. 



C. Participation number 

The qualitative features of the fringe-visibility dynamics, given some initial preparation, are determined by its 
participation number defined as: 



M = 



n -1 



= number of participating levels in the LDOS 



(26) 



In the case of a TwinFock preparation, the Wigner function is spread all over the equator of the spherical phase-space 
and thus overlaps with all the states in the sea up to the scparatrix level (the equator intersects with all sea trajectories 
in Fig. [2] but with no island trajectory). Therefore we expect M to be of order N, with classical (N independent) 
prefactor that reflects the relative size of the sea: 



M 



ClassicalPrcfactor x N, [TwinFock preparation] 



(27) 



In the case of a coherent preparation, the Wigner function is a minimal wavepacket that has width cr„ = (iV/2)^/^. 
the Fock states \n) were the eigenstates of the Hamiltonian, as they are in the Fock regime, we would get 



If 



M = (27riV) 



1/2 



[Coherent preparation, n basis] 



(28) 



for all three coherent preparations. However in the Josephson regime the eigenstates are \E^), and therefore the 
differences between the LDOS of the three coherent states Zero, Pi, and Edge, are set by the ratio between cr„ and the 
width of the separatrix {An = NK/U). The ratio cr„/An equals the dimensionlcss semiclassical parameter (u/NY/"^. 
If this ratio is larger than unity the distinction between the Zero, the Pi and the Edge preparations is blurred, and we 
expect to have the same participation number. Indeed plotting the participation number of these three initial states 
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(FiglU, we see that at the strong interaction Umit M (3/2)iV-'^/^. This is roughly half compared with Eg. ((28)) . and 
reflects the odd-even selection rule that removes half of the overlaps (the semiclassical states do not have a well-defined 
parity with respect to site substitution, whereas the actual quantum eigenstates are constructed from their odd and 
even superpositions. The participating constituents are those that have the same parity as that of the preparation). 
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FIG. 4: (Color online) The participation number M as determined from the LDOS for N = 100 (o), 500 (□), and 1000 (o) 
particles. The left panel contains the Zero (lower set in blue) and Pi (upper set in red) preparations, while the Edge preparation 
is presented in the right panel. Note the different vertical scale. In the crudest approximation we expect in the Edge case 
M ~ A^^'^^, while in the Pi case M < iV^/^ as long as {u/N) < 1 



For weaker interaction, when the semiclassical parameter {u/Ny^'^ is smaller than unity the different nature of the 
Zero, the Pi and the Edge preparations expresses itself. Now we have to account for both the width of LDOS line 
shape and the mean level spacing. This leads to the following results: 

[Zero preparation] (29) 

[Pi preparation] (30) 

[Edge preparation] (31) 

The striking point here is that the Pi preparation resembles the Zero preparation, rather than its sister separatrix Edge 
state. This seems at first site in contradiction with semiclassical intuition: one would naively expect that wavepackets 
that have the same energy and reside in the same phase space region (separatrix) would behave similarly. This is not 
the case, as we see here, and later in the dynamical analysis. The Pi state is actually closer to the Zero preparation, 
as both have a small participation value. The resemblance of the Pi and Edge preparations is only detectable in the 
formal limit N — )■ oo. In other words, because of the N dependence of AI it is "easier" to approach the "classical 
limit" in the case of an Edge preparation. 
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V. COHERENCE DYNAMICS 



A. Classical, semiclassical, and quantum dynamics of the Bloch vector 

Mean field dynamics that is generated by using the GPE, describes the classical evolution of a point in phase space, 
which is a single trajectory. The semi-classical theory describes the classical evolution of a distribution in phase 
space. The Quantum theory adds to this picture recurrences and fluctuations which are due to the discreteness of the 
energy spectrum. The dimer system is intergrable, and therefore the WKB method provides a very good basis for the 
analysis. Fig|5] illustrates the agreement between the quantum evolution of the Wigner function, starting from the 
TwinFock state [l^, and the semiclassical evolution of a corresponding distribution. It should be emphasized that 
in the Wigner- Wyle formalism any operator A is presented by the phase-space function Aw(i^), and the calculation 
of an expectation value can be done in a classical-like formulation: 

(a) = trace[p M = j ^Pw(^^)^w(^^) (32) 
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FIG. 5: (Color online) The evolving quantum state of = 40 bosons with u = 5 for TwinFock (n — 0) preparation. The units 
are such that K = 1 and the time is t = 4. On the left - the Wigner function of the evolved quantum state. On the right - the 
corresponding classical evolution. 

The singlc-particlc density matrix of the two-mode system can be presented using the Bloch vector, 

S={J)/{N/2) = {S^,Sy,S,) , (33) 

Whereas previous work |^ has focused on the mean- field dynamics of the occupation difference {N/2) (Sz), here we 
study the single-particle coherence manifested in the one body purity, 

OneBodyPurity = (1/2) [l + {S,)^ + {Sy)^ + {S^f] , (34) 

and the transverse component of the Bloch vector, the fringe visibility. 

In Figj6]we plot examples for the evolution of the Bloch vector, and observe significant differences in the dynamical 
behavior of the Zero, Pi, and Edge preparations. Our objective is to understand how the dynamics depends on the 
dimensional parameters: the "classical" parameter u and the "quantum" parameter N. 




t (scaled) t (scaled) t (scaled) 



FIG. 6: (Color online) The variation of Sx{t) with time for A''=40 particles with it = 5, for Zero (left). Pi (middle), and Edge 
(right) preparations. Note the different vertical scale. The dashed-dotted lines are based on semiclassical simulation. 

In accordance with the opening paragraph of this section, we see in Figj6] that in the semiclassical simulation the 
fluctuations always die after a transient. This should be contrasted with both the classical (single trajectory) behavior, 
and the quantum behavior. In the latter quantum case the wavepacket is a superposition of M > 1 eigenstates, and 
consequently there are persistent fluctuations that depend on the "quantum" parameter N. In what follows, we 
will quantitatively analyze the characteristic features of the quantum dynamics resulting from the three coherent 
preparations, including the frequency of oscillation, its mean long-time value, and its RMS amplitude, and compare 
the analytic predictions to the numerical results at various values of the characteristic parameter u/N . 
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B. Characteristic frequencies 



The typical frequency of the fluctuations is the simplest characteristic that differentiate the three panels of FigJB] 
The numerically-obtained frequency w^^^ as a function of the interaction parameter u/N is displayed in Fig .[71 In 
the classical picture uj^^^ should be related to the Josephson frequency wj at the bottom of the sea, while quantum 
mechanically it reflects the level spacings of the participating levels. A straightforward analysis leads to the followings 
estimates: 

[Zero] (36) 
[Pi] (37) 
[Edge] (38) 
[u > N] (39) 

The first three expressions apply for (u/N) < 1, where the differences between the preparations is distinct. The 
last expression refers to the regime {u/N) > 1 where the differences are blurred and the three preparations become 
equivalent equatorial states (as clearly evident from Fig.[7l). Note that due to the mirror symmetry of the Zero 
preparation the expected frequency should approach 2wj, while for the Pi preparation it is bound from below by 2u)^. 
Both frequencies are indicated in Figl7]by dashed lines. 
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FIG. 7: (Color online) The mean frequency of the Sx{t) oscillations versus u/N for for A*" = 100 (o), 500 (□), and 1000 (o) 
particles. The preparations are (upper to to lower sets of data points): Zero (blue), Edge (magenta), and Pi (red). The 
weaker-interaction theoretical predictions (|36|l - (|38| l. which are doubled due to mirror symmetry, are represented by blue, red, 
and magenta dashed lines, while the strong interaction prediction H39p is represented by a black dash-double-dotted line. 



C. Long time average 

Next we examine the long time average, plotted in Fig|5] as a function of the characteristic parameter. Unlike the 
fluctuations around it, this average value does not reflect the quantization of energy and therefore a purely semiclassical 
analysis is adequate. The naive expectation might be that coherence would be diminished due to phase spreading 
(aka "phase diffusion"). This is indeed the case in the fock regime, leading to {Sx)oo ~ 0. However, the situation 
is rather more complicated in the Josephson regime, where {Sx)oo is determined by u/N. The semi-classical phase 
space picture allows to calculate the phase distribution P{ip) that pertains to the long time ergodic-like distribution 
(see for example FigE]). This distribution is determined by the LDOS. Then we use the integral 

J cos{ip)P{ip)d^ , (40) 

to evaluate the residual coherence. This procedure results in the following predictions: 

1/3 [TwinFock] (41) 

51 w exp[-(u/iV)] [Zero] (42) 

^« _l„4/log[^(u/iV)] [Pi] (43) 
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Thus, the coherence of the Zero preparation is robustly maintained as long as u/N < 1, corresponding to the (f = 
phase locking of the two condensates due to the weak coupling. By contrast, the Pi and Edge coherence is far more 
fragile throughout the Josephson regime [l^l- Note that the phase locking of the Pi state would take place only if 
it were the Rabi-regime (no separatrix). As evident from the Twin-Fock self- induced coherence (|^T|) . one-particle 
coherence should not necessarily be lost due to interactions, but could actually be built This course of events is 
somewhat similar to the coherent relaxation of a system to its ground state at low temperatures: the ground state 
has higher purity compared with the initial preparation. 



i< p a a op OD tEnga. ' Q:) 




u/N u/N 



FIG. 8: (Color online) Left: The long-time average of S^t) versus u/N for N = 100 (o), 500 (□), and 1000 (o) particles. The 
preparations are (upper to lower sets of data points): Zero (blue), TwinFock (black), Edge (magenta), and Pi (red). The 
symbols are used for the quantum results and the dashed lines are the semiclassical prediction for fourty particles. Note that 
the scaling holds only in the Josephson regime 1 <^ u <^ , and therefore, for a given u/N range, becomes better for large A'^. 
Right: The long time RMS of Sx{t) for the three coherent preparations (lower to upper sets): Zero (blue). Edge (magenta), 
and Pi (red). In the inset, the RMS of Sx{t) for Edge (A) and Pi (V) preparations is plotted versus A*' while u = 4 is fixed. 



D. RMS of the fluctuations 



Finally, we turn to discuss the RMS of the fluctuations, which constitute a fingerprint of energy quantization. 
General reasoning implies that the classical fluctuations are suppressed by factor M: 



RMS [(A) J = 



I C4uj)duj 



-,1/2 



(44) 



In this formula Cci(a;) is the classical power spectrum of an crgodic trajectory. What wc want to highlight is the 
quantum N dependence. It is important to clarify that in the semi-classical limit M — !■ oo, and therefore the 
fluctuations are suppressed. We emphasize again that the quantum behavior is intermediate between the coherence- 
preserving classical (single trajectory) dynamics and the strong coherence-attenuation of the semi- classical (infinite 
M) dynamics. 

The dependence of M on N is remarkably different for the various preparations. In the case of the TwinFock 
preparation M ^ N and therefore the RMS is inversely proportional to 7V^/^. This should be contrasted with the 
case of the Pi and the Edge preparations: 



RMS [S^it)] 
RMS [S^{t)] 
RMS [S^{t)] 



7V-1/2 

iV-i/4 

(log(/V))-l/2 



TwinFock] 

Edge] 

Pil 



(45) 
(46) 
(47) 



We further note that in the Pi case the leading dependence of the participation number is on the classical parameter 
(M ^ u^/^), unlike the case of the Edge preparation where the leading dependence in on the quantum parameter 
{M ^ 7V^/^). The results of the RMS analysis are presented in Fig[51 The implied 7V^/^ scaling based on Eq. (Hil) ) is 
confirmed. The dashed lines in the inset are power-law fits that nicely agree with the predictions of Eq. (|45p . 
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VI. SUMMARY 

Using semiclassical machinery, we have analyzed the temporal fluctuation of the single-particle coherence, and of 
the fringe-visibility, in the Bose Josephson model. While recent experiments in the Josephson regime have essentially 
focused on mean-field population dynamics @, with coherent preparations at the spectral extremes (self trapping 
versus Josephson oscillations), here we highlight intricate effects that can be found by studying the coherence-dynamics 
in the intermediate separatrix regime. We predict significant differences in the transverse relaxation of seemingly 
similar coherent initial states, differing by the initial relative-phase or by their location along the separatrix, as well 
as the interaction-induced phase-locking of two initially separated BECs due to the combined effect of interaction and 
coupling. The semiclassical WKB quantization facilitates the calculation of the LDOS for the pertinent preparations, 
and thus the estimation of the number of participating eigenstates. This allows for a detailed quantitative analysis of 
the time evolution of the fringe visibility function. The obtained analytic expressions are found to be in a very good 
agreement with the results of numerical calculations. 
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